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In this paper we derive quark model results for scattering amplitudes 



meson contains a heavy quark. This "1313" system is an attractive theo- 
retical laboratory for the study of the nuclear force between color singlets; 
the hadronic system is relatively simple, and there are lattice gauge flie- 



s' ■ ory (LGT) results for VeB(r) which may be compared to phenomenological 

models. We find that the quark model potential (after lattice smearing) 
has qualitative similarities to the LGT potential in the two B*B* channels 
in which direct comparison is possible, although there is evidence of a dif- 
ference in length scales. The quark model prediction of equal magnitude 
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but opposite sign for 1=0 and 1=1 potentials also appears similar to LGT 
results at intermediate r. There may however be a discrepancy between 
the LGT and quark model 1=1 BB potentials. A numerical study of the 
two-meson Schrodinger equations in the (bq)(bq) and {cq){cq) sectors with 
the quark model potentials finds a single BB "molecule", in the 1=0 BB* 
sector. Binding in other channels might occur if the quark model forces are 
augmented by pion exchange. 

I. INTRODUCTION 

The origin of the residual strong force between hadrons is a complicated problem. 
Several distinct scattering mechanisms have been suggested as important contributors 
to interhadron forces, and it may be difficult to distinguish these experimentally. As 
an example, models of the NN force have been proposed which include t-channel meson 
exchanges, short-range quark-gluon interactions, intermediate s-channel excitation of A 
baryons, and various other effects. Comparisons with NN data alone may not determine 
the relative importance of these mechanisms, since one might find an unphysical param- 
eter set that happens to describe the data well with a particular scattering mechanism, 
especially if there are many free parameters. 

This complication is illustrated by a "confusion theorem" which notes that the two 
mechanisms most often assumed in models of the NN force, t-channel meson exchange 
and quark interchange, can easily be misidentified since they correspond to identical flavor 
flow. Both scattering mechanisms are of course present in nature, and the problem is to 
determine their relative importance as a function of separation. One can see that they 
are physically distinct because they represent scattering through intermediate states in 
different sectors of Hilbert space, one additional qq pair for meson exchange versus no 
extra pairs for quark interchange. 

Lattice gauge theory provides an attractive opportunity to isolate the contributions 



of the various mechanisms that have been proposed for residual interhadron forces. By 
taking the limit of a very heavy "6" quark and introducing sources for B = bq mesons, 
one can study the energy of a meson pair as a function of separation |l]-§[. (We follow 
and use B generically to refer both to a pseudoscalar B and a vector B*; technically bq is 
an anti-B meson, but the results for scattering amplitudes and potentials are identical.) 
A lattice B meson has a fixed heavy-quark coordinate, and in the BB system one can 
use this to determine the energy Ebb{t) of the BB pair as a function of center-of-mass 
separation. The difference between this energy and that of two isolated B mesons pro- 
vides a natural definition of the Vs#(r) interhadron potential. By changing the initial 
and final coordinates of the light-quark Green functions one can in effect vary the light 
quark flavor, and thereby determine the identical BB (actually 1=1) and distinguishable 
BB (1=0) potentials. Changing the meson source angular quantum numbers allows one to 
infer separate BB, BB* and B*B* potentials, providing that the associated multichannel 
mixing ambiguities can be resolved. One may also investigate the importance of different 
scattering mechanisms by evaluating potentials associated with different quark lines di- 
agrams, such as direct versus quark interchange. Finally, in the more difficult full-QCD 
simulations one can test the importance of additional qq pairs in hadronic forces. Clearly, 
many questions which are of great importance to model builders may be answered by this 
application of lattice gauge theory. 

In this paper we evaluate the various BB potentials in the context of the nonrelativistic 
quark model, for comparison with existing and future LGT results. At present, configu- 
ration mixing in LGT constrains the direct comparison to two B*B* channels, but there 
is already evidence of qualitative agreement. Statistically more accurate LGT results and 
separation of the various B and B* spin and isospin channels should allow very interesting 
comparisons with the various BB potentials we derive here. 

The paper is organised as follows: Section II introduces the Coulomb plus linear 
quark model and the technique used to evaluate hadron-hadron T-matrix elements, and 



carries out the detailed evaluation with SHO wavefunctions. Section III gives the general 
relation between the T-matrix and equivalent local potentials, and uses this to derive the 
BB potentials. Section IV discusses the details of these BB quark model potentials and 
compares these to LGT results, and studies the possible formation of bound states using 
the BB Schrodinger equation. Finally, Section V gives a summary and conclusions. 

II. BB T-MATRIX FROM THE QUARK-GLUON INTERACTION 

A. Method and previous applications 

The technique we use to determine quark model Vgg(r) potentials is to evaluate the 
lowest (Born) order T-matrix element of the interquark Hamiltonian between two-meson 
scattering states, which is then Fourier transformed to give an equivalent low-energy 
potential. The interaction assumed is the OGE color Coulomb and spin-spin interaction 
and linear scalar confinement. The effective interquark hamiltonian for this interaction is 
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where the sum runs over all pairs (i,j) of valence quarks and antiquarks that are in 
different initial hadrons. (Pairs of quarks in the same hadron contribute to hadron energies 
rather than to scattering; the partition of H into Hq and Hi in this formalism is well- 
known in atomic physics, and is discussed elsewhere in the hadronic context |7|||.) The 
color generator in Hi is as usual T a = X a /2 for quarks and T a = —X aT /2 for antiquarks. 
After a single interaction of this Hi between a constituent pair in different initial hadrons, 
quark line interchange is required to give an overlap with the color-singlet final meson 
states. For meson-meson scattering this gives four diagrams, which are shown in Fig.l. 

This model (incorporating only the spin-spin OGE term, which is dominant in light 
hadrons) gives an excellent description of S-wave meson-meson scattering in channels 
without valence annihilation, specifically 1=2 tttt pj and 1=3/2 Kit ||, with 2 and 3 pa- 
rameters respectively (a s /m%, Psho', m,/m s ). These successful results are impressive in 



that the parameter values are already well known from light meson spectroscopy, and the 
optimum values found in fitting the scattering data alone are consistent. This suggests 
that, at least for PsPs scattering at moderate energies, Born-order quark-gluon diagrams 
with external meson wavefunctions describe the dominant scattering mechanism. KN 
S-wave scattering at low momenta is also excellently described by this model. (A good 
simultaneous fit to higher- momentum S-wave KN scattering however requires a some- 
what reduced nucleon wavefunction length scale [ 10| ; this may be due to short- distance 



correlations in the nucleon's three-quark wavefunction, which is not included in our simple 
Gaussian forms.) In all these successfully modelled reactions there is of course no one- 
pion-exchange term, since a three-pseudoscalar vertex is not allowed. There are also no 
s-channel resonance contributions; these specific reactions were studied precisely because 
they do not have the complication of valence qq annihilation. Studies of the NN inter- 
action in the quark model JTT|, using both perturbative and nonperturbative techniques, 



have found a large short-range repulsive NN core interaction due to this OGE interaction, 
and similarly conclude that the dominant core interaction at short distances arises from 
the OGE spin-spin hyperfine term. 

B. Evaluation of BB scattering amplitudes 

In this paper we evaluate the contribution of all three terms in Eq.(l) to the BB T- 
matrix elements {T/j} and potentials {Vrb(t)}. These will be presented with separate 
flavor, color, spin and space factors, so the BB case can easily be generalized to other BB 
spin channels. Although the spin-spin hyperfine term was found to be dominant in light 
pseudoscalar-pseudoscalar S-wave scattering, in BB we also expect the color Coulomb 
interaction to be important, since it must dominate at short distances. We will derive the 
{Tfi} for general quark masses, with fa the light (antiquark) mass and m the heavy "6 
quark" mass. 

To evaluate the meson-meson potential we first calculate the matrix element of the 



quark Hamiltonian Eq.(l) between two-meson initial and final states. Conservation of 
three-momentum implies that this matrix element is of the form 

(CD\H Z \AB) = -±-T fi 5(A+ B - C - D) . (2) 

(In previous scattering calculations we instead gave results for a Hamiltonian matrix 
element hfi, which is trivially related to the Born-order T-matrix element by hfi = 
Tfi/(27i) 3 .) In our earlier discussion of 7T7t scattering we distinguished the four scat- 
tering diagrams according to which pair of constituents interacted; these are "capturei" 
(CI), "capture 2 " (C2), "transferx" (Tl) and "transfer/ (T2); see Fig.l. The hadron- 
hadron Tfi matrix element for each diagram can conveniently be written as an overlap 
integral of the meson wavef unctions times the underlying quark-level T-matrix element. 
In the quark-quark Tfi (Fig. 2) the initial and final constituent momenta are a,b — * 
a',b'. It is useful to write the quark-quark Tfi in terms of the linear combinations q, p± and 
p 2 , defined by q = a' — a = b — b', pi = (a + a')/ 2 and p 2 = (b + b')/2. For the specific case 
of one gluon exchange, the complete quark-quark Tfi to second-order in three-momenta 
(suppressing the color factor) is 
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This follows from taking the matrix element the one-gluon-exchange effective hamiltonian 
^A^j" between an initial and final quark pair, and using the definition Eq.(2) of Tfi. 
We have displayed the 7070 terms and the 7^ terms separately in this Tfi] the 7j7« terms 
are proportional to \/mim 2 . This result is valid for both quarks and antiquarks; only 
the color factor distinguishes them. For completeness we give the corresponding T-matrix 
element due to linear scalar confinement, which is 
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In the four overlap integrals that result from taking the two-meson matrix elements of 
these quark T-matrices (corresponding to the four independent scattering diagrams) we 
find that f>2 is constrained to equal ±pi plus a diagram-dependent shift. These overlap 
integrals are explicitly (introducing A = (m — in) / \m + m) , and using p = pi as an 
integration variable) 

T ( £ l \AB ^CD) = 

ffd 3 qd 3 p <5>* c (2p + q- (1 + X)C) <5>* D {2p- q- 2A - (1 - X)C) 
T fi (q,p, -P + C) <5> A (2p - q - (1 + X)A) <5> B (2p- q- (1 - X)A- 2(7) , (5) 

TJf 2) (AB -► CD) = 

JJd 3 qd 3 p %(-2p + q + 2A - (1 + A)(7) $^(-2p - <f- (1 - A)(7) 

7>(<f,p, -p - C) $ A (-2p + q+ (1 - X)A) $ B (-2p + q+ (1 + X)A - 2(7) , (6) 

t£ 1] {AB^CD) = 

ffd 3 qd 3 p <5>* c {2p + q- (1 + X)C) ®* D {2p - q- 2A-{1- X)C) 
T fi (q,p,p-1-C) $ A (2p - q- (1 + X)A) $ B (2p + g - (1 - X)A - 2(7) , (7) 

t£ 2) (AB^CD) = 

JJd 3 qd 3 p %(-2p + q + 2A- (1 + X)C) ®* D {-2p - q- (1 - X)C) 
T fl (q,p,p-A + C) $ A (-2p + q + (1 - X)A) $ B (-2p - <f + (1 + X)A- 2(7) . (8) 

With standard quark model SHO wavefunctions (given in Appendix A) each overlap 
integral above becomes the quark Tfi times a shifted Gaussian. The overlap integrals are 
then (also assuming elastic scattering in the CM frame, so \A\ = \C\) 
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d 3 qd 3 pexp\ - — (p - p ) 2 \ exp | - — (q - q f \ T fi (q,p, -p - C) , (10) 
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The shifts p~o an d <fo are diagram dependent, and are 

' q/A + (A + C)/2, CI 



(q-qoy\T fi (q,p,p-A + C] 
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(13) 



(14) 



q/4 + {A-C)/2, C2 

(A + C)/2, Tl 

1(1- C) /2, T2, 

'2(-A + AC)/3, ClandC2 

;i + A)(-l+C)/2, Tl 

-(l-X)(A + C)/2, T2. 
Note that the CI and C2 integrals over p must be carried out before the q integral, since 

the po shift depends explicitly on q in this case. 

These results are for a general quark Tfi(q, pi, P2). In this paper we consider the special 
case of ^-independent quark interactions, corresponding to pure V(r,ij) quark potentials 
in coordinate space. This simplification is appropriate for the color Coulomb, linear 
scalar confinement and spin-spin contact hyperfine interactions treated here; the spin- 
spin hyperfine term simply has an additional multiplicative spin factor for each diagram. 
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This assumption is not valid for the spin-orbit and tensor interactions, which have explicit 
p dependence; these will be treated in subsequent work. Given a quark T-matrix of the 
form 



T f i(q,Pi,P2) =T fi (q) 
we can further simplify the SHO overlap integrals above. This gives 
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where % for each diagram is given by Eq.(14). 



C. Explicit meson-meson T-matrix elements 



(19) 



We will now evaluate these overlap integrals with the quark Tfi due to color Coulomb, 
spin-spin hyperfine and scalar confinement interactions, in Eq.(3) and Eq.(4) and trans- 
form these into equivalent low-energy Vbb potentials. The specific quark interactions we 

use are (with color and spin factors removed) 

' Ana s /q 2 , color Coulomb 

Tfi(q) = I —(8na s /3mimj), spin-spin hyperfine (20) 

k 6irb/q 4 , linear confinement. 



9 



1. Spin- spin Hyperfine Contribution 



The spin-spin hyperfine contribution is derived using the overlap integrals above and 
the color and spin matrix elements given in our previous discussion of 1=2 tttc scattering . 
The results are presented as the meson-meson T- matrix element 7+ iagr (AB — > CD) = 
(signature) • (flavor factor) ■ (color factor) • (spin factor) • [space overlap]. We also define 
a frequently occurring combination 



IT 2 = (1 - \f{A + C) 2 + (1 + \) 2 {A- C) 2 . 



The results are 



7> (C1) = (-1) • (1) • (-4/9) • (-3A 

T fl {T1) = (-1) • (1) • (+4/9) • (3/8) 
T fl {T2) = (-!)■(!)■ (+4/9) -(3/8) 
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(A-C) 



(21) 



(22) 
(23) 
(24) 

(25) 



2. Color Coulomb Contribution 



The four color Coulomb overlap integrals can be evaluated similarly using the quark 
color Coulomb Tfi in Eq.(20), which gives the results 



7> (C1) = (-1) ■ (1) • (-4/9) • (1/2) 
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(27) 

(28) 
(29) 
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TABLE I. Physically Allowed BB States. 
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gular Quantum Numbt 
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Stot = o 


Stot = 1 


Stot = 2 


BB 
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even L 


- 


- 









oddL 


- 


- 


BB* 




1 


- 


allL 


- 









- 


allL 


- 


B*B* 




1 


even L 


oddL 


even L 









oddL 


even L 


oddL 



3. Scalar Confinement Contribution 



Finally, with the linear scalar confinement Tfi we find 



7v(ci) 
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(30) 
(31) 

(32) 

(33) 



D. T-matrix elements for physical BB states 

Since the BB and B*B* systems have identical mesons there are constraints on the 
physically allowed states; these are summarized in Table I. The physical BB scattering 
amplitudes are diagonal in isospin, since we have assumed equal light quark masses. 
To extract these isospin-diagonal amplitudes we evaluate the T-matrix element between 
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BB pairs with definite isospin, for example \B B ) for 1=1. With our phases the B 
(b) and B (6) meson isodoublets are {\B°), \B~}} = {-\bd),\bu)} and {\B+),\B°)} = 
{ — \ub), — \db)}, analogous to the kaon system. Since \B~B~) = \(bu)(bu)), this implies 
identical antiquarks as well as quarks. As noted in Ref. J7|, in this case there is a second set 
of "symmetrizing" quark line diagrams, with quark lines exchanged rather than antiquark 
lines. These have the effect of interchanging the final mesons C and D; when added to 
the antiquark exchange diagrams of the previous section this gives a Bose-symmetric 
scattering amplitude, satisfying the even-L constraint in Table I. For BB (or BB) the 
complete 1=1 BB elastic scattering amplitude is then 

t bb (i=i) = Tf .( AB _^ CD ^ + T fi (AB -> DC) , (34) 

where Tfi(AB —>■ CD) is the sum of Eqs. (22-33) of the previous section. Similarly for 1=0 
BB we find a second, symmetrizing diagram, but with an opposite sign; 

t bb (1=0) = _ Tf .( AB ^ CD) + T fl (AB -► DC) . (35) 

This gives a spatially antisymmetric scattering amplitude. Thus 1=1 BB is allowed only 
even L and 1=0 BB is allowed only odd L. Another consequence of the relative signs in 
T^ ~ ' above is the relation between Born-order BB potentials in systems that differ 
only in total isospin, 

V& = °\r) = -Vg =1 \r). (36) 



III. V B b potentials 

A. Potentials from the T-matrix: general formalism 

A2^2 T-matrix can be represented as an equivalent Born-order potential operator 
V op .{x\ — a>2, Vi, V2), between pointlike particles 0]. The definition of this potential 
operator is 
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S(A + B- C- D) T fi (A,B,C,D) = 

d 3 x l d 3 x 2 e- i{d -^^ 2) V op [x x - x 2 , Vi, V 2 ) e +l{A - £l+§ -* 2 > . (37) 



(2tt) 3 

To evaluate this potential operator for a given T-matrix one can write the meson-meson 
T fi (AB ->■ CD) as a function of the variables Q = (C - A), P 1 = (A + C)/2 and 
P2 = {B + D)/2. A power series expansion in the {Pi} variables is then performed, 

T fi (AB - CD) = T^\Q) + T^°\Q) P u + T^\Q) P 2l + T^\Q) P u P 2j + ... . (38) 

The {Pi} in the T-matrix expansion are replaced by left- and right-gradients in the equiva- 
lent potential operator defined implicitly by Eq.(37). This procedure gives a local potential 
operator that reproduces the specified scattering amplitude Tfi at Born order. One may 
confirm that this approach reproduces the full 0(v 2 /c 2 ) Breit-Fermi Hamiltonian when 
applied to the one photon exchange e~e~ Tfi, expanded to 0(P 2 ). 

The leading term T^°\Q) in the Pi expansion is a function of Q only, and Fourier 
transforms into a local (static limit) potential that is a function of X\ — x 2 — r only. 
In the cases we consider here T^(Q) is a function of \Q\ only, which leads to a local 
potential that is a function of r only. The relation between Tf t (Q) and V(f) is 

V(f) = j±y 3 jd 3 Q T<°>(<?) e*5> . (39) 

In this paper we obtain a local Vbb(^) potential by Fourier transforming the BB — >■ BB 
T^(Q), which we obtain from the full meson-meson Tfi by changing to the variables 
{Q, Pi} and setting P 1 = P 2 = 0. 

B. Vbb : individual contributions 

1. Color Coulomb 

As an illustration we shall evaluate the 1=1 BB potential due to the color Coulomb 
interaction in the infinite m limit (A = 1), following which we will simply quote the re- 
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maining results. The 1=1 B B Tfi matrix element we find for A = 1 with this interaction 

is 



T 



BB (7=1) 



fi 



2 3 7TO S 



% lFl (l/2, 3/2; ^-^-) - lFl (l/2, 3/2; £z£)l , - | 



3 1 / 



6[3 2 



2(3 2 



-(A-C) 2 /2f3 2 



+ (C^ -C) 



(40) 



This is the sum of Eqs. (26-29) for A = 1, symmetrized as in Eq.(34). In this case there is 
an obvious partition into "direct meson" and "crossed meson" scattering contributions; 
the direct contributions have a forward-peaked Gaussian in Q 2 = (A — C) 2 . Since the 
direct Tfi is a function only of Q 2 no expansion in Pj is required, and we obtain the 
potential simply by Fourier transforming; 



V&-"(r) 



£.« J *Q e*+*»» {^ lFl (l/2,3A J) -.F.(l/2,3/2;g 

Evaluation of these integrals is discussed in App.B; the result is 



(41) 



V, 



(J=i), 



_BB 



color Coulomb 



2a, 
9r 



(2/vr) 1/2 /3r-4Erf(/3r/2) 



-/3 2 r 2 /2 



(42) 



The three contributions in the square brackets are from Tl, T2 and (C1+C2) respectively. 
The small-r behavior has an obvious interpretation; for r = r^ much less than the wave- 
function length scale f3~ l , the Born-order heavy quark-quark interaction term Tl must 
approach the bare color Coulomb result — 2a s /9r. The remaining quark- ant iquark and 
antiquark-antiquark terms retain mean constituent separations of 0(/9 _1 ) as r — > and 
so have nonsingular limits. 
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2. Spin- Spin Hyperfine 



In the limit of large quark mass only the T2 diagram has a nonzero Tfi with this 
interaction, which for 1=1 BB is given by Eq.(25) with A = 1. Since this is a simple 
Gaussian, the Vbb{t) resulting from Eq.(39) is also a Gaussian, 



V, 



(J=i), 



BB 



2 1 ' 2 a s [3 3 



spin— spin hyperfine 



9tt 1 /2 



-f3 2 r 2 /2 



m z 



(43) 



3. Linear Confinement 

The T-matrix elements of the linear confining interaction for 1=1 BB are given by 
Eqs. (30-33). In the A = 1 large quark mass limit the forward-peaked part of Tfi equals 
Anb 



Tfi 



3 1 / 2 iFxC-1/2, 3/2, ^) + ^(-1/2, 3/2, |£) + 1 



2 /2/3 2 



• (44) 



direct ' 

Evaluation of the Fourier transform of this Tfi requires an integral which is discussed in 
App.B. The result for Vbb(j') is 
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(45) 



We have again grouped terms according to diagram. The first square bracket gives the Tl 
(quark-quark) term, the second is T2 (antiquark-antiquark) , and the third is the rather 
complicated C1+C2 quark-antiquark term. 

As with the Coulomb overlap integrals we could have anticipated some properties 
of this potential. First, at small r the interaction of two heavy quarks approaches the 
bare Vftb( r w>) times a color and spin factor of 2/9 (instead of the usual qq color-singlet 
coefficient 4/3). Thus the Tl potential approaches ((2/9)/(4/3)) -br = br/6 for r « f3~ l . 
The antiquark-antiquark (T2) and quark-antiquark (C1+C2) potentials again approach 
finite limits at small r, and give a contact potential of 
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The largest individual diagram contribution at contact is the positive T2 (antiquark- 
antiquark) term; the mean antiquark-antiquark separation is larger than quark-antiquark, 
which gives a larger linear-potential matrix element. However there are two contributing 
quark-antiquark diagrams, CI and C2, which give equal contributions; their sum is larger 
than T2 and opposite in sign, so at contact we find a net attraction. At larger r the sign 
of this interaction is reversed. 



C. VaMr) final results 
The full 1=1 BB potential is given by 
v£-"(r) = -£ { 1 + PA) 1 ' 2 Pr ~ 4 Erf(/Jr/2) } e"**" + |^ ^ e^' 2 

+ |. { Pr ,-*+!* + ^ «-^/» - (>r + -I) Erf(/3r/2) e -^/ 2 _ _i_ e ^V/4 j , 

(47) 

which is the sum of the color Coulomb, OGE spin-spin and linear confinement contribu- 
tions. 

The potentials for the remaining BB potentials can be obtained similarly. In all cases 
we find that for 1=0 there is a simple relative flavor factor which changes the overall sign 
of Vbb-i as i n Eq.(36). The various B*B* potentials can be determined from Vbb above by 
changing spin overlap matrix elements, which are given in Table II. For example, to convert 
the Vgg~ (r) potential in Eq.(fJ?D to Vg*g*' tot ~ (r) one multiplies the color Coulomb and 
linear contributions by (+l)/(+l/2), and the remainder, the spin-spin hyperfine term 
(ex a s /m 2 ), by (+l/4)/(+3/8). 

The BB* potentials require more careful treatment. Just as we found in BB, the 
BB* T-matrix has forward- and backward-peaked contributions, but they are no longer 
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TABLE II. hot = 1 BB Spin Matrix Elements (from Table I of 



System 



Operator 



mesons 


•SW 


o = i 





— <->a • Of, 


BB 





+1/2 




+3/8 


BB* -> BB* 


1 


+1/2 




+3/8 


BB* -> £*£ 


1 


+1/2 




-1/8 


B*5* 


2 


+1 




+1/4 




1 







+1/2 







-1/2 




+5/8 



identical in magnitude; this was required for BB by Bose symmetry at the meson level. It 
is again useful to associate these with a "direct" BB* — > BB* potential (from the forward- 
peaked contributions to the T-matrix) and a "crossed" BB* — ► B*B potential from the 
backward-peaked contribution, in which there is a B <-+ B* transition at each "crossed- 
V" interaction. At Born order in S-wave scattering the direct- and crossed-potentials can 
just be added to give a total effective BB* potential. This total S-wave BB* potential 
has twice the BB Coulomb and linear potential and +2/3 of the BB spin-spin potential, 
which makes it identical to the S to t — 2 B*B* potential. The spin matrix elements for 
the direct and crossed BB* contributions are given in Table II. 



IV. DISCUSSION 



A. Numerical results for V BB (r) 



We show the total V BB (r) of Eq.(47) and the three individual contributions in Fig.3. 
The parameters employed are a s = 0.5, b = 0.18 GeV 2 and rh = 0.33 GeV, which 
were used by Scora and Isgur in their recent HQET discussion of B meson semileptonic 
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decays |13| . They quote several values of the variational best-fit SHO j3 for B mesons, 



specifically (3 = 0.43 GeV (5), 0.40 GeV (B*) and 0.35 GeV (IP B 3 mesons); we adopt 
an intermediate value of 0.40 GeV. 

The total 1=1 BB quark model potential is evidently strongly attractive at small r, 
passes through a node at r « 0.28 fm, and is weakly repulsive at larger r. The Coulomb, 
spin-spin hyperfine and linear confinement contributions to Vg^~ (r) for r ss 0.5-1 fm are 
all repulsive and are comparable in magnitude. 

The short-range attraction is dominantly due to the color Coulomb attraction; for 
r « (3~ 1 the bound-state wavefunctions are irrelevant, and we see an unscreened color 
Coulomb potential between the heavy quarks, with a color-spin factor of 2/9. This gives 
an attractive short distance potential H&( r ) — ~~ 2a s /9r. At small r this quark-quark 
interaction diagram Tl is dominant; at larger r the other color Coulomb diagrams and 
bound state screening become important, and the Coulomb contribution crosses over to 
a weak repulsion at r « 0.36 fm. The Coulomb contribution is +8 MeV at 0.5 fm, and 
by 1 fm it has fallen to +2 MeV. 

Of course at sufficiently small r the OGE Born approximation will be inaccurate, and 
the bbqq system will deform to minimize the dominant small-r color Coulomb interaction. 
In 1=1 the most attractive channel has bb in a color 3; this should give a stronger color 
Coulomb force than our Born result, and with these deformed wavefunctions our Born- 
order relation V^ I=0 ^ = — V^ I=1 ^ will not be accurate. 

The contribution of the linear confining interaction to Vgg (r) is not large because 
there are approximate cancellations between the four diagrams and (unlike Coulomb) 
there is no regime in r in which one diagram dominates. We find that the linear contri- 
bution to V^b i s attractive at short distances, with a contact value of about —50 MeV, 
and crosses over to a weak repulsion at r « 0.38 fm. At 0.5 fm the linear confining term 
contributes +7 MeV and at 1 fm it is +6 MeV. 

In light-quark hadrons such as 1=2 tttt and NN (the core potentials) one finds that 
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the color spin-spin hyperfine term makes the dominant contribution to the hadron-hadron 
interaction. Here we instead find that at moderate r the hyperfine, linear and Coulomb 
potentials make comparable contributions. The smaller hyperfine contribution to the 
BB system follows from the absence of both capture diagrams and one transfer diagram; 
these vanish due to the 1/m.jmj prefactor in Hj Eq.(l). The spin-spin capture diagrams 
in particular made the largest contribution to the 1=2 tctt interaction. We find that the 
hyperfine contribution to V BB ~ (r) is repulsive (as in 1=2 tcti), but is much smaller; the 
contribution to V^b (r = 0) is +26 MeV, which falls to +16 MeV at r = 0.5 fm and 
+3 MeV at r = 1 fm. 

B. Comparison to LGT BB potentials 

Several references have discussed the determination of BB potentials using LGT tech- 
niques [0-0. The most detailed study to date is by the UKQCD collaboration ||. This 
work treats the b quark as a static, spinless source, so there are four potentials, labelled by 
the light antiquarks' total isospin and spin, (J t *f ,S t l Jt ) = (0,0), (0, 1), (1,0) and (1, 1); 
these are shown in Figs. 6-9 of Ref. ||. Both (0,0) and (1,1) show strong short-range 
attraction. (1,1) appears consistent with weak repulsion beyond contact (the first lattice 
point is at r « 0.18 fm). The (0, 0) potential shows a clear rise to a (probable) zero near 
0.3 fm, and some evidence for weak repulsion at larger r. The (1, 0) and (0, 1) potentials 
are small (ca. ± 50 MeV) and are not yet well characterized, although (0, 1) shows some 
intermediate-range attraction, and both potentials show evidence of repulsion at contact. 

Comparison of our quark model BB potentials to these LGT results is unfortunately 
nontrivial except for S to t = 2 B*B*, due to the spin degree of freedom. The lattice static- 
quark limit has degenerate B and B* mesons, so the lowest-energy configuration for given 
S^ot will not be the external source basis state (such as \BB), as in our quark model 
calculation) but instead will be the linear superposition of \BB), \BB*) and \B*B*) that 
gives the lowest expected energy |TJj . 
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A direct comparison does appear possible for the UKQCD S t l 9 t = 1 potentials, which 
correspond to an S to t — 2 B*B* meson pair. These have no S-wave degenerate BB or B*B 
mixing states, and should therefore be similar in physical meaning to our S to t = 2 B*B* 
quark model potentials, provided that the tensor coupling to L=2 BB is unimportant. 
In Figs. 4, 5 we compare the (1,1) and (0,1) LGT potentials to our 1=1,0 S tot = 2 B*B* 
potentials. Clearly there is qualitative similarity, although the quark model potentials 
appear to have a larger length scale. 

A more realistic comparison is possible if we apply an estimated lattice resolution effect 
"smearing" to our quark model potential. Lattice resolution can have a dramatic effect on 
some aspects of the potential, especially near r = where there is little Jacobean weight. 
For example it will regularize a continuum 1/r term, so that the LGT BB potentials 
approach finite values at contact, as noted by Stewart and Koniuk |Q. To model lattice 
resolution effects we introduce a Gaussian-averaged quark model potential, defined by 

V(r)= fdh'—^— exp{-(r-r')7a 2 }VV) . (48) 

We choose the smearing length a to be the lattice resolution of 0.18 fm estimated by 
UKQCD for their LGT results ||. The resulting V potentials are shown as dashed lines 
in Figs. 4, 5. Except for what appears to be a difference in the length scale, these are 
qualitatively similar to the LGT potential. In future we should ideally compare with 
LGT potentials from simulations that have a finer spatial resolution. 

The isospin dependence of the quark model BB potentials is a very characteristic 
feature. The 1=0 Born-order quark model potentials are equal in magnitude but opposite 
in sign to the 1=1 potentials. This result may be supported by LGT at intermediate r 
(compare the LGT (0,0) with -(1,0) and (0,1) with -(1,1) in [g); (0,1) and (1,1) are also 
shown in our Figs. 4, 5. At contact however the LGT 1=0 and 1=1 results differ greatly in 
magnitude; since 1=0 is odd-L it may be difficult to extract the small-r 1=0 BB potential, 
and in any case we expect the Born result to be inaccurate at small r, because the strong 
color Coulomb term will dominate. 
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C. Bound States 

Bound meson pairs, known as "molecules" , most easily form in channels in which the 
pair can exist in S-wave. From Table I the S-wave BB channels are 1=1 BB, 1=0,1 BB*, 
and (I, Stot) = (1, 0), (0, 1) and (1, 2) B*B*. We have searched for possible bound states in 
these BB and T>T> systems by numerically integrating the two-meson Schrodinger equation 
with the generalizations of Eq.p7|) to the different spin channels. 

With our standard parameter set (in Sec. IV. A) and assuming B and B* masses of 
5.288 GeV and 5.325 GeV, we find that only one channel has sufficient attraction to form 
a bound state; this is 1=0 BB*, which has a deuteron-like repulsive core and intermediate- 
range attraction. We find a binding energy of just —5.5 MeV with our parameters, typical 
of nuclear binding energies. 

The most attractive of the 1=1 attractive-core channels are BB* and Stot — 2 B*B*, 
which have identical potentials. The attraction however is not strong enough to overcome 
the intermediate-range repulsion. As we increase a s we do find that these systems bind, 
but at an unphysical a s ~ 1.0, about twice the usual quark-model value. 

In all the T>T> channels the smaller reduced mass makes binding more difficult, and we 
find no bound states. 

One pion exchange forces are often suggested as an important component of the meson- 
meson interaction, and have been discussed in general by Tornqvist [15] and quantitatively 



by Ericson and Karl [jig] . Ericson and Karl find that one pion exchange is not attractive 
enough to bind mesons lighter than BB, but that the 1=0 Stot = 2 (odd-L) B*B* channel 
will bind from this interaction alone. The S-wave BB channels with attractive one pion 



exchange forces are 1=1 Stot = B*B* QT5 1 and 1=0 BB* jL7|. Since we expect both 



one pion exchange and quark-gluon forces to be present in nature, one might study the 
combined effect of the one pion exchange potential and the analogues of our Eq.(47) in 
a search for bound states in other channels that are more accessible to experiment than 
BB. 
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V. CONCLUSIONS 

We have calculated the T-matrix and low energy equivalent potentials between pairs 
of heavy- light mesons, the "£>£>" system, in the context of the nonrelativistic quark model. 
The assumed scattering mechanism is a single interaction of the standard quark model 
Hamiltonian, with OGE color Coulomb and spin-spin terms and linear confinement. The 
parameters used were taken from previous studies of meson spectroscopy and HQET 
matrix elements. This model of the hadron-hadron T-matrix is known from previous work 
to give a good description of experiment in the analogous light pseudoscalar channels 1=2 
mi and 1=3/2 Ktx. We carry out the overlap integrals of this interaction with standard 
SHO external bq meson wavefunctions in closed form, and so obtain analytic results for 
Vg B ~" ' (r) in the various allowed channels. These are compared to recent LGT results 
from the UKQCD collaboration in the channels where this is possible, which are 1=0,1, 
S to t — 2 B*B*. We find results similar to these UKQCD potentials after lattice smearing, 
but with a somewhat larger length scale. Our 1=1 BB potential however is attractive at 
small r, which appears inconsistent with UKQCD results. 

We find that our quark model potentials are sufficiently attractive to support a bound 
state in only one channel, 1=0 BB*, which has a deuteron-like potential. With pion 
exchange added, other channels may have sufficient attraction to support bound states. 

In future work a more detailed comparison with LGT potentials, especially below 
0.2 fm, will be important as a test of the forces assumed in the quark model calculation. 
It would also be very interesting to generate LGT BB potentials for large but finite 
quark mass, so the meson spins could be specified uniquely. One could then compare 
the LGT and phenomenological BB potentials in all channels unambiguously. Finally, 
one may extract the spin-dependent (spin-orbit, spin-spin, tensor and so forth) equivalent 
BB meson channels using using similar techniques, and we anticipate that a comparison 
with LGT results for these potentials might also be interesting as a test of the nature of 
spin-dependent "nuclear" forces. 
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APPENDIX A: MESON WAVEFUNCTIONS 

In this appendix we present the explicit meson wavefunctions used in the paper. A 
single meson state is given by 

\A(A,S,S,)) = 

E -^=5 c5 J2(S,S z \l/2,s z ;l/2,7s z ) f fd 3 a d 3 a 5(A- a - E) $ A (3 rel ) | <&,£-,) > ( A1 ) 

c,g=l V 3 s z ,s z J J 

where the relative momentum variable is a re i = (ma — ma)/((m + m)/2). The full 
momentum-space wavefunction is 

& A (Aa,a) = 5(A- a-a) & A (a re i) (A2) 

We normalize this state to 

(A(A,S,S Z )\A(A',S',S' Z )) = 6(A-A')6ss>8s,s>, , (A3) 

and the individual quark and antiquark states are similarly normalized as 

(q(a,s,s z )\q(tf,q' z )) = 5(a-a') 8 SzS > z . (A4) 

This implies a relative-momentum wavefunction normalization of 

r d 3 (a rel /2)\$ A (a rel )\ 2 = l . (A5) 
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The full and relative spatial wavefunctions are related by 

^A(x cm , r) = , 2 , 3/2 Mr ) (A6) 

with x cm = (mxq + rhXq)/(m + rh) as usual. These are related to the momentum-space 
wavefunctions by 

$ A (A,a,E) = -ip JJd 3 x cm d 3 r e -^+S). gcm -ia rRl .r/2 ^( fcm y) . (A7) 

The relative spatial wavefunction i^a{t) is similarly related to the relative momentum 
wavefunction $^4 by 

Mr) = T^Ll Jd\a rel /2) e +i ^ 2 > r ® A (a rel ) , (A8) 

where r — x q — Xq. 

The ground state SHO quark model wavefunction which we use in the potential cal- 
culations discussed in the text is 

®A(a rd ) = $ (a rel ) = 3/403/2 ex P { ~ drel V 8 / 3 *} • (A9) 



APPENDIX B: OVERLAP INTEGRALS 

In deriving BE potentials as Fourier transforms of the scattering amplitudes we en- 
countered shifted-Gaussian overlap integrals of the form 

I a ,c = Jd 3 Q <r**lP«** {Ffa c; cS^) . (Bl) 

To evaluate integrals of this type it is useful to introduce an integral representation of the 
confluent hypergeometric function. For c > a > we use 



1 F 1 (a,c;x) = — -^ f ' dt e xt f 1 (1 - tf 

1 [a) lie — a) Jo 



y-*- 1 (B2) 

l\a) l\c — a) Jo 
which leads to 
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r(o) T(c — a) \c J c\ 



■n \ 3 / 2 r a 



dyy a - 1 {y + lfl 2 - c {k-yy- a -\ 



-py 



(B3) 



where k = (co/ci — 1) x and p = /3 2 r 2 /4co. For the color Coulomb interaction we have 
c = 3/2 and c = a + 1, so the integral becomes 

/i A 3/ 2 = /rf 3 Qe-^ 2 /^ + ,QV iFi(1/2! 3 /2;Ci ^ ) = Z^^ f kp ds s-^ ( 



P 2 ' 24c 1 / 2 P 1/2 Jo 



2/32 



7T 2 (3' 



(cqCx)- 1 / 2 e~ /32r2/4co Erf (c 2 (3r) 



(B4) 



where C2 = c (co/ci — 1) 1 ^ 2 /2. The special case Co = ci 



is 



A/2,3/2 



^^ ^/3 2 r 2 /4c 



c =ci CoT 

The linear confinement overlap integrals lead to the case a = —1/2 and c 
required integral is 



(B5) 
3/2. The 



Z-i/2,3/2 = d*Q e-«*"f 



lFl (-l/2,3/2; Cl ^) 



Q 2 

75 2 



7T 2 tf 3 r 1/2 



2 c 2 



/3 2 r 2 /4c 



|( Co / Cl -l)V2( w+ i.) Erf(w) _ 
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+ 



7T 



3/2^3 



2(c - Cl ) 3 / 2 



3 -/3 2 r 2 /4(c - Cl ) 



(B6) 



and the special case c = C\ is 

-7-1/2,3/2 



2 /34 
7T P __ _-/? 2 r 2 /4c 



Co=Cl 4c 
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Fig.l. The four meson- meson scattering diagrams. 
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Fig. 2. Momentum definitions in the quark-quark T-matrix element. 
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Fig. 3. The V^ B potential, showing individual contributions. 




Fig. 4. Comparison of the Vg~ B * tot ~ quark model potential 
(solid is calculated, dashed is smeared by a = 0.18 fm) with the 
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LGT potential of Ref. 
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Fig. 5. Comparison of the V B , B * tot quark model potential 



(solid is calculated, dashed is smeared by a = 0.18 fm) with the 
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(flight— 0, Slight — 1) 
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LGT potential of Ref. 
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